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Abstract 

We show that Ma’s scotogenic model, which is arguably one of the simplest settings 
containing a Dark Matter candidate and generating a naturally suppressed active 
neutrino mass at 1-loop level, suffers from a potentially severe hierarchy-type prob¬ 
lem. In case the right-handed neutrinos involved have sufficiently large masses, they 
can via loop effects drive the mass parameter of the inert scalar contained in the 
model towards negative values. This behaviour leads to a breaking of the Z 2 parity 
symmetry built into the model which is paramount to keeping the setting consistent 
at low energies - without it the model would lose its Dark Matter candidate and 
the neutrino mass would not be naturally suppressed. Thus, if the breaking occurs 
at a sufficiently low scale, it could potentially spoil the success of the whole model. 
Trying to avoid this consistency problem leads to a new constraint on the model 
parameter space which has not yet been described in the literature. 
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1 Introduction 


The smallness of neutrino masses and the identity of the Dark Matter (DM) are among 
the greatest puzzles of modern particle physics, as is the question of how to extend our 
Standard Model (SM) and how to probe the extensions. The most attractive theories 
beyond the SM are those which can address several of the known open problems. This 
is often the case for models generating a light neutrino mass only at loop level (i.e., 
radiatively) - cf. e.g. Ref. Q to see how such settings are constrained at both, low and 
high energies. Depending on the particle content, there exist models which generate an 
active neutrino mass at 1-loop [^, 2-loop [^|^, or 3-loop level, but probably the 
simplest extension compatible with all data is Ma’s scotogenic model (^. 

The scotogenic model just adds three right-handed (RH) neutrinos Nji and a scalar 
doublet Tj to the SM, all of which are charged under an additional Z 2 parity symmetry]^ 
This symmetry is crucial for the model to work; without it, neutrino masses would already 
be generated at tree level and none of the possible DM candidates of the model could 
be stable. However, if the Z 2 is intact, the scotogenic model cannot only account for 
phenomenologically valid neutrino masses jS 10 and potentially for DM it 

can also lead to a variety of phenomena in low-energy experiments such as lepton flavour 
and/or number violation (^[9,14-17 or in high-energy collider searches 18 21 , as well 


as to new aspects for neutrino model building [Io|[22|[2^ . 

A few years ago, the hrst study of the renormalisation group running of the scotogenic 
model has appeared in Ref. 24 . It has been shown that several parameters exhibit strong 


running effects. This is known in similar settings: Ref. 25 showed that an inert scalar 


can trigger electroweak symmetry breaking (EWSB) at 1-loop level even if not present at 
tree level - a fact that remains true if the setting is extended 26 or discretised (2^. More 


generally it is well-known that scalar mass parameters are sensitive to large scales 28 


which is also true for models with two Higgs doublets 29 


In this paper we apply a similar logic, i.e., we investigate corrections to the mass 
parameter of the inert scalar. The heavy RH neutrinos can drive the squared mass 
parameter of the inert scalar towards negative values, thereby giving a non-zero vacuum 
expectation value (VEV) {rj) and by that breaking the crucial Z 2 parity symmetry, which 
is why we call this observation the parity problem of the scotogenic model. This could be 
disastrous given that, if the parity was not a conserved global quantum number, not only 
would the model lose the stability of its DM candidates but also a neutrino mass would 
be generated at tree level at phenomenologically unacceptably large values. 

On the other hand, given that the breaking happens at a high scale, the question is 


^Due to the Z 2 symmetry restricting its interactions, the new scalar doublet is usually called inert. 
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justified why we should bother for low-energy phenomenology. The answer is two-fold: 
hrst, broken discrete symmetries could potentially lead to problems in cosmology. If the 
Z 2 symmetry was e.g. intact at a high scale but then at some point broken by the large RH 


neutrino masses, we would potentially be in danger of creating unwanted domain walls 30 


which could modify the history of the Universe in an undesirable way. At the same time, 
a VEV of the inert scalar would lead to large tree level masses for active neutrinos. 
This would cause them to immediately freeze-out and possibly lead to an intermediate 
matter-dominated phase of the Universe, which would again alter the expansion history. 
However, one could argue that the domain wall problem can be cured 31 35 and that 
active neutrinos may simply re-thermalise as soon as the Z 2 is intact again. In addition, 
the symmetry will be restored at sufficiently high temperatures 36 . Thus, given that 


all these processes happen before big bang nucleosynthesis, there may not even be an 
observable remnant. This is a perfectly justihed viewpoint, but even if one disregards 
the points above it turns out that the breaking scale of the Z 2 can be as low as a few 
TeV, as we will show in Sec. 4T If that is the case, DM production can be signihcantly 
modihed in the scotogenic model and points in the parameter space which are thought to 
be consistent and to lead to successful DM production could in fact be in trouble due to 
the parity problem. Ultimately, the message is that one has to be careful and check for 
any given parameter point whether a potential Z 2 breaking leads to a problem, or not. 

As a side note, it is interesting to mention that obviously any issues associated with 
a too large RH neutrino mass scale could be avoided if all RH neutrino masses were 
sufficiently small. From this point of view, one could even argue that to some extend the 
scotogenic model would have a preference for light sterile neutrino masses. Since they 
would then also be lighter than all neutral scalar components, they would automatically 
be the actual DM candidates of the model. Such settings are known to work very well (see 
e.g. Refs. 37 50 for suitable production mechanisms), and they have also been discussed 
from a phenomenological point of view in the context of the scotogenic model [8,13. 
While it is not easy to hnd a suitable mechanism to motivate light sterile neutrino masses 
in the scotogenic model 51 , the above arguments could be interpreted as such a scenario 
with light sterile neutrinos being in fact quite natural. 

This paper is organised as follows. After a brief overview of the model in Sec. 
we discuss in Sec. [^ the general possibilities for its possible vacuum conhgurations and 
illustrate the approximate constraints arising from avoiding a violation of the Z 2 parity 
in the scotogenic model, which yields a simple but accurate formula. A more advanced 
numerical analysis of a few detailed examples is presented in Sec. We discuss some 
of the aforementioned implications for cosmology in Sec. [^ before we conclude in Sec. 
Technical details can be found in the Appendix. 
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2 Model Overview & Constraints 


The scotogenic model is one of the most minimal frameworks combining a naturally 
small neutrino mass at 1-loop level with several DM candidate particles. The particle 
content is that of the SM, supplemented by (typically) three RH neutrinos {i = 1, 2, 3) 
and a second scalar doublet rj with SM quantum numbers identical to those of the Higgs. 
The crucial addendum is an additional Z 2 (parity) symmetry, under which all SM particles 
are neutral whereas the new fields are odd. It is this symmetry which simultaneously leads 
to the light neutrino mass being generated at 1 -loop level only and to the stability of the 
potential DM candidates. 

Compared to the SM, several qualitatively new terms appear in the Lagrangian: 

• The RH neutrinos get a direct Majorana mass term + h.c., which leads 

to masses Mi 2,3 upon diagonalisation. 

• A neutrino Yukawa coupling ^Yukawa D + h.c. {fj = ia 2 r]*) involving the 

new scalar and the RH neutrinos in addition to the SM lepton doublets It is 
crucial to observe that this term does not lead to a tree level neutrino mass, as long 
as the Z 2 is unbroken and thus prevents the held r] from obtaining a VEV. 


• The scalar potential involving the SM Higgs 0 and the rj held is given by: 
V = - 1 - 

-h As (0V) + -^4 (hV) (0M + y [(^7^)^ + ^-c.] . 


( 1 ) 


In this expression both mass parameters mf 2 must be real, as need to be the cou¬ 
plings Ai_ 2 , 3 , 4 ; As can be chosen real and positive by absorbing its phase into rj. 


Note that it is the combination of the Majorana mass term, the new Yukawa coupling, 
and the As-term in Eq. ([^ which violates lepton number. If any of those coefficients 
was zero, a global U{1) lepton number could be consistently dehned and the symmetry 
of the Lagrangian would be increased. Thus, by virtue of’t Hooft naturalness 52 , the 


renormalisation group equations (RGEs) for those quantities will only allow for changes 
proportional to the couplings themselves, so that they remain small everywhere if they 
are small at any energy scale, cf. the Appendix. 

The scalar potential ([^ needs to yield EWSB without compromising the Z 2 parity. 

This suggests the parametrisation 0 = ^0,n -|- and rj = {'r)^-,rpY'^ where only the 
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physical fields are written down explicitly. Splitting into real and imaginary parts, 
^ {Vr + gives rise to the following physical scalar masses 


53 


2\iv\ 

(2a) 

ml + ^'^A 3 , 

(2b) 

^2 + (A3 + A4 -|- As), 

(2c) 

^2 + (A3 + A4 — As). 

(2d) 


Note that the parameters in the scalar potential are subject to a number of theoretical 
constraints originating from the requirement that the scalar potential be bounded from 


below 54-56 


Ai > 0 , A2 > 0 , A3 > —\/AiA2 , A3 + A4 — IA5I > —A1A2, 


(3) 


which also affects the valid mass spectra resulting from Eqs. (2a) to (2d). In addition, we 


can impose the condition that we remain in a perturbative regime, that is Ai, 2 , 3 , 4,5 
In our numerical analysis, we demand that the above conditions remain valid for the 1- 
loop corrected running couplings up to a given high energy scale. The 1-loop RGEs can 
be found in the Appendix. 


3 Vacuum Structure New Constraints 

The scalar sector of the scotogenic model is a particular realisation of a two Higgs doublet 
model (THDM) with a Z 2 symmetry imposed on some fields. In a general THDM both 
scalar doublets (called and 02 here to ease the comparison) may acquire a VEV and we 
find the tree level vacuum of the theory by minimising Eq. ([^. Replacing the fields by their 


VEVs ^ 1 ^ 2 , which can be chosen real if both electric charge and CP are conserved 54,59 
leads to two equations: 


vi {m\ + \iv\ -I- Xvl) = 0 and V2 {m\ + + A't'i) = 0, 

where we have abbreviated A = A3 -|- A4 -|- A5. 


(4) 


^Note that there are some ambiguities in the definition of perturbativity 57 58 . For definiteness, we 
have chosen the criterion Ai, 2 , 3 , 4,5 < Itt in our numerical computations. 


4 

















These equations allow four qualitatively different sets of VEVs0 

@ uH = 0, 

^ 2 _ Azm^-Am^ o _ Aim^-Amf 

— A2-A1A2 ’ ^2 — A2-A1A2 • 

While (T) and (2) respect the Z2, (3) and (T) break it spontaneously. In our numerical 
analysis, we have investigated whether (?) or (T) are stable minima of the potential. Oth¬ 
erwise the Z 2 would be broken in any case since the potential must have a minimum and 
only Z 2 breaking vacua are left. 

Computing the Hessians we hnd that we have a stable, Z 2 symmetric vacuum if 


ml > 


0 if > 0 ( 1 ), 

if < 0 @. 


(5) 


The latter equation is equivalent to the condition m|, > 0, i.e. the held rjji develops a VEV 
if its mass square - an eigenvalue of the Hessian - becomes negative. Had we allowed for a 
relative phase between vi and U 2 , we would hnd that Im 02 (i-e. ///) may develop a VEV. 
The condition to exclude this is ml > 0. Similarly, to avoid breaking electric charge, 
we need m\ > 0. If we ignore a possible instability of the Z 2 symmetric vacua, we may 
expand the theory around the wrong vacuum. Expanding around the correct vacuum, 
however, would alter the phenomenological predictions of the model and compromise the 
Z 2 symmetry. 

Avoiding vacuum instability, of course, does not exclude that the Z 2 symmetric minima 
are not the global minimum of the potential. We can investigate this by plugging the 
solutions (?) ^ @ back into the scalar potential, Eq. ([^, which yields 61 : 


% = K = 0,^2 = 0) = 0, 


I vt = 






V(^^V{vt = 0,vl = -j^ 


4 

nil 


nq 

' 2 A 2 


(6a) 

(6b) 

(6c) 


^( 2 ) and @ each have two real solutions, which are physically equivalent due to global symmetries, 
while (D gives rise to four solutions, which reduce to two physically inequivalent solutions 60 . 


5 









V®^V vt = 


2 X2'mi — Xm\ 2 Aim^ — \m\ 


A2-A1A2 


-,^2 = 


A 2 - A1A2 


\ 2 TrLl — 2\mfm2 + Aim2 
2(A2-AiA2) 


( 6 d) 


Any of the above solutions (T) - (7) is the global minimum if Vj 2 ^0 and if it has the 
lowest vacuum energy, Vq < V® for all fc 7 ^ h Note that, however, there is some freedom 
in choosing a minimum which is not the global one if the decay time is larger than the 
age of the Universe, see e.g. 60 . Nevertheless, we can use the above equations to check 


whether the breaking of Z 2 is the only constraint, or if there could be further constraints, 
which originate from the minima (T) and ( 2 ) not being the global one. 


For illustrative purposes, we now show analytically that the Z 2 symmetry can be 
broken spontaneously due to radiative corrections. We do so by directly calculating the 
breaking scale. To simplify the equations, we only consider one generation of fermions 
and assume m\{^) > 0 for all scales /r, i.e. we study under which conditions only 
becomes negative. For sufficiently small quartic scalar couplings we can ignore their 
contributions to the RGBs altogether and consider two simple but illustrative limiting 
cases, where all quantities under consideration are assumed to be real. In addition, we 
shall assume that the Majorana mass M does not run at all, which is justihed numerically, 
cf. Ref. 


24 


Case 1: Large neutrino Yukawa coupling — h{fi) ^ 

In this limiting case, neglecting any gauge and non-neutrino Yukawa couplings, we can 


approximate the coupled RGBs (A-2b) and (A-5b) for h and respectively, as follows: 

5 


Vh ~ -h ^, 
2 ’ 


— 2h?‘m\ — 


(7) 


where V = and we have suppressed explicit scale dependence for 

brevity. These differential equations can be solved exactly and yield a condition for the 
mass square to become negative, i.e. for symmetry breaking to occur at some point f = f*. 


=0 




-(■-Sf) 


5/2 


( 8 ) 


This yields in the linear approximation [where mKO) -C 2M^]: 

dvr^m^O) 


f ~ 


M2h2(0) ■ 


(9) 
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Case 2: Small neutrino Yukawa coupling — h{^) (y'j(/i) 

In this case the RGBs involve the gange conplings Eq. (A-1): 


~ {gl + 3 ^ 2 ) h, Vml ~ {g^ + 3gl) mj - 

The solutions are simple to find and we get 

47r^m2(0) 


( 10 ) 


i* — 


M2/l2(0) 


( 11 ) 


Note that the conditions for both cases agree if m^O) <C 2M^, which is just the case we 
are interested in. 

One might object that a Z 2 broken above some large scale is irrelevant as long as 
the low energy phenomenology is unaffected, in which case we should only be concerned 
whether G is large enough for this to be true. On the other hand, one may argue that a 
broken symmetry in the UV is undesirable from an aesthetical point of view, and we can 
impose that the crossing to negative values of shall not occur up to some high scale. 
In fact, both requirements are not so different, because if we require the Z 2 symmetry 
to be intact at least to the TeV scale Mxev = 10^ GeV (in order not to disturb DM 
productions) or to, say, the scale of “grand unification” (GUT) Mgut = 10^® GeV (where 
new physics may modify the situation), the corresponding bounds for po = Mz are: 


7712(0) > 


In (103 GeV/Mz, 10^6 GeV/Mz) 




h(0)M ^ (0.246, 0.905) Mh(0). (12) 


In both cases, as long as the tree level scalar mass parameter 7772 ( 0 ) is not at least of 
the size of the heaviest Majorana mass, radiative Z 2 breaking could potentially occur. 
However, since we at the same time do not want 7772 ( 0 ) to be too large in order not to pull 
the Higgs mass parameter to too large values and in the worst case even spoil EWSB, cf. 


Eq. (A-5a), some tension between both requirements is generated. 


In the upcoming section we will discuss both the appearance and the impact of radia¬ 
tive Z 2 breaking. 


4 Numerical Analysis 

In this section, we will show by a numerical example how non-trivial constraints on the 
model parameter space can be obtained if the requirement of keeping the Z 2 intact is 
imposed on some range of energies. We thereby use the full system of RGEs for all three 
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generations. The goal of this section is to illustrate that keeping the Z 2 intact up to a 
certain energy scale leads to non-trivial bounds. 

Depending on the setting the Z 2 breaking may be problematic if it, e.g., occurs at 
scales where DM production could be modihed. However, we would like to stress that, 
due to the size of the parameter space, it is not possible to make any global statements. 
Rather, whenever a parameter point is studied, one would need to check on a case-by-case 
basis how/if phenomenology is affected by the resulting constraint. 

In general, the model parameters “beyond the SM” are M^, m2, A2, A3, A4, A5}, 
some of which are constrained by low- or high-energy data. Nevertheless, we treat these 
numbers as free parameters, to keep the computational time limited and since e.g. a 
particular leptonic mixing does not affect our considerations. We £x A 5 = 10“® at the 
input scale /i = Mz, whose value is kept small by the corresponding RGE (A-4e). The 
conditions ([^ require that we use large enough values for A 2 , such that it is positive 
for all energy scales. We hnd that, at p = Mz, A 2 = 0.1 gives good results while 
the numerics are under control. For the neutrino Yukawa couplings hij, we input a 
bimaximally mixed setup of 0(0.1) at the GUT scale, which is known to potentially 


yield phenomenologically valid leptonic mixing at low energies 24 


Moreover, the heavy helds are integrated out for renormalisation scales below their 


mass thresholds 63 and the SM input values have been chosen according to 64 


The connection between the different parameters is intricate: at the input scale, we 
must choose m\ large enough for it not to be driven to negative values by the corrections 
from the Majorana masses, cf. Eq. (A-5b). At the same time, a large m\ will generate large 
positive corrections to ml [provided that 2 A 3 -|- A 4 > 0, see Eq. (A-5a)], which could spoil 
EWSB. This is where the actual tension resides: the running mass parameter m\ must 
lie within the range of the electroweak scale and at the same time it must not be too small 
to avoid Z 2 symmetry breaking. There is one more player in the game: the RH neutrino 
masses drive the tension by their appearance in the RGE for m^, Eq. (A-5b). If chosen 
too large (say, around 10 TeV, for hij ~ 0.1), demanding an unbroken Z 2 symmetry at 
all scales requires ml{Mz) to be larger than allowed by the obligation to achieve EWSB. 

As mentioned in Sec. the inert doublet VEV may not only lead to breaking of the 
Z 2 , but it could also cause electric charge and/or GP violation. Note however that, since 
A 5 is small for all renormalisation scales, we will have nearly degenerate GP-even and -odd 
neutral scalar masses, mu ~ m/, such that GP violation is inseparable from the breaking 
of Z 2 . To unambiguously identify such scenarios, we have made use of the equivalence 
of a squared mass becoming negative to symmetry breaking (see Sec. |^. This has been 
investigated by replacing the Lagrangian parameters in the mass relations (2b) to (2d) 
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Figure 1: Left: Z 2 -breaking scale A as a function of the masses m±^Fi for R-H masses 
(Ml, M 2 , M 3 ) I= (900,1500,5000) GeV. Right: Explicit running of tur {< m±) 
for one example set of parameters with (Mi, M 2 )|^^^^^^ = (900,1500) GeV (dashed grey 
lines). In both plots, Yukawa couplings < 1 are chosen (the crucial ones are given). 


by the running couplings and the tree level VEV by a running VEV (cf., e.g.. Ref. |^): 


2 2/ ^ 


(13) 


This can be done because, in the broken phase, we only need counterterms that are 
invariant under the broken symmetry group 66 . Thus, we can obtain the counterterms 


in the broken phase from those in the symmetric phase. We can use this fact to construct 
running quantities in one phase from running quantities in the other 67 . 


4.1 Breaking scale 


Let us first try to get a feeling for the numbers involved, i.e., for the scale where the Z 2 
breaking occurs as a function of the physical scalar masses]^ 

The natural question is at which scale the breaking could occur for realistic particle 
masses. The crucial observation is that even RH neutrino masses in the TeV range could 
be perfectly sufficient to break the Z 2 parity at a relatively low scale, as Eqs. ([^ and ( [II| 
suggest. To illustrate this statement, we present in the left panel of Fig. [T]the breaking 
scale as a function of the scalar masses m± and uir. We have neglected all SM fermions 
but the top quark to be able to determine the breaking scale with higher precision. The 


^We could have used the scalar mass parameters in the Lagrangian, but we consider physical masses 
to be more illustrative. On top of that, one can use Eqs. (2b) to (2dI to convert one into the other. 


9 




















RH neutrino masses are {Mi, M 2 , M^) = (900,1500,5000) GeV, input at Mqut, with 
Yukawa couplings of 0(0.1), which are all not overly large and more or less in the ranges 
where many particle physicists would “naturally” expect them. However, glancing at the 
hgure, it is visible that even for these values the breaking scale can be as low as a few 
TeV. The breaking scale grows with both m± and mn, as to be expected from Eqs. (2b) 
and (2c), since both of them grow with m 2 but are related as = (A4 + As)^^. 

Ultimately, the value of the breaking scale is determined by a tug-of-war between the 
scalar mass parameter m 2 and the maximum of all {Mihij}, cf. Eq. (12). 

To see the exact influence of the combination of Yukawa coupling and mass, we can 
look at the right panel of Fig. [T| where the explicit evolution of m^ is shown as a function 
of the renormalisation scale, for three different exemplary masses M 3 . Indeed, as soon as 
M 3 becomes dynamic, it pulls mji towards smaller values, until it is zero. Beyond this 
scale m\ is negative, i.e. the “physical” mass m^ becomes imaginary. However, given 
that this occurs in a setting with m\ < 0 and m^ < m±, the scale at which mR = 0 is to 
be identihed with the Z 2 breaking scale, so that in fact the assumed held conhguration is 
no longer expanded around a minimum of the potential for larger renormalisation scales. 
Expanding around an actual minimum would of course yield real and positive masses. 

An important point to make at the end this section is to understand that it is not 
immediately visible whether or not the breaking emerges at sufficiently low scales and 
destroys the validity of the model. Instead, one has to look at the details of the parameters 
under consideration, and if in doubt apply further analyses that allow to clarify the 
situation, as we shall outline in Sec. However, as can be seen from the right panel of 
Fig-B it is well possible that e.g. for an inert scalar with a mass of 262 GeV, which is the 
DM candidate of that setting if it is the lightest electrically neutral and stable particle 
involved, the symmetry breaking can already occur at about 1 TeV, which is very close to 
the phase decisive for DM production in the early Universe. In such a case, the running of 
the parameters involved cannot simply be neglected when computing the DM production. 
This observation is also true if the effect of Z 2 breaking is related to thermal effects, which 
may signihcantly alter the situation 68 -72 , as we will explicitly discuss in Sec. 


4.2 A humble way out 

For some studies it may not be easy to show whether or not a too low breaking scale 
has a negative impact, e.g., if no suitable tools to compute DM production are available]^ 
In other contexts one may only be interested in a rough picture of the situation, where 

^While DM production can be handled easily in simple cases, in particular settings with co¬ 
annihilations or non-thermal production are more involved. 
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the decisive point is to avoid a potential symmetry breaking altogether. In such cases, a 
convenient way out would be to require the breaking scale to be at such high energies that 
it does not affect low-energy physics. For example, one could ask for which parameter 
values the Z2 breaking happens only beyond Mqut, where new physics would potentially 
modify the situation. In the following, we give an example for how this would affect the 
parameter space available. 

To illustrate the parity problem, we generated a random set of 10® input values within: 


OTeV < m 2 {Mz) < 3TeV, —1 < < 1, 


(14) 


for masses {Mi, M 2 , M 3 ) = (900,1500,2000) GeV as input at Mgut- We then computed 
the running effects for each point chosen, this time taking into account all SM fermions. 
First, we solve the RGBs for the running couplings of the scalar sector (Ai, 2 , 3 , 4 , 5 ) and check 
the scalar consistency criteria, Eq. (|^, for all energy scales between Mz and MgutI^ If 
these are not violated for any scale below Mgut, we solve the scalar mass parameter 


RGBs (A-5a) and (A-5b). If we do not find EWSB below 1 TeV, i.e. if m^(p) > 0, we 
reject the input values. Otherwise we distinguish two sub-cases depending on the sign of 
m\{jj), see Eq. ([^. In case we find that all vacua that respect the Z 2 parity are unstable, 
we also reject the input values. By adhering to this approach, we can avoid rejecting 
input values because of a broken Z 2 symmetry although some other criterion fails as well, 
i.e., we determine exactly those points which would not be rejected if one disregarded the 
bound arising from radiative symmetry breaking. 

The results of our parameter scan are shown in Fig. where the allowed and rejected 
input values at the electroweak scale are marked]^ All data points violating vacuum 
stability or perturbativity have been excluded from the plot. The colour code is as follows: 
Black dots mark the points which fulfil all the constraints. All other points are excluded 
for various reasons. Yellow dots are excluded due to failing to produce EWSB below 
1 TeV, while red/green dots are excluded only because they lead to radiative Z 2 breaking, 
where rn^j^ < 0 is signified by red dots and m/_ < 0 by green ones. 

For a large range of input values for A 3 and A 4 , the criterion that we must encounter 
EWSB below 1 TeV translates into an upper bound on m 2 {Mz), which is a function of 
A 3 and A 4 . Only if 2 A 3 -(- A 4 < 0, EWSB can be maintained below 1 TeV while allowing 
for a large m 2 . This manifests itself in the columns of valid (black) data points ranging 
up to very high input values of m 2 , which are visible in both panels of Fig. 


®This ensures that we do not have to consider further particle thresholds and all SM fermions can be 
safely neglected. This even holds for the top quark given that its mass is so close to the lower input scale. 

^This time, we have chosen to display the Lagrangian parameters as opposed to Fig. since the 
resulting bound is illustrated in a better way. 
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m 2 [GeV] 


Figure 2: Parameter scan for the scotogenic model. Black dots are valid input parameter 
values, yellow dots violate EWSB, red < 0) and green (m^ < 0) dots are excluded 

due to breaking of Z 2 . The shaded grey area represents the “lower bound” on m 2 . 

On the other hand, one can see that including the criterion of an unbroken Z 2 symme¬ 
try is essentially equivalent to a lower bound on m 2 {Mz) - in our example about 550 GeV, 
indicated by the grey area in the plots. As can be seen in the left Fig. we have found 
red and green points below this bound, which are only excluded by the requirement of 
the Z 2 symmetry to be unbroken at all scales. Even above this bound it is visible that, 
for certain choices of A 3 ^ 4 , there are some red and green dots which would be considered 
unproblematic if parity breaking was disregarded. 

For the physical scalar masses in the minimal subtraction (MS) scheme, we obtain 
values in the ranges 


554.2 GeV < m±(/x = m±) < 2780.6 GeV, 


558.5 GeV < mR/i^n = rriR/i) < 2781.6 GeV, 


(15) 


conhrming our expectation of a lower bound on the masses. Generalising our estimate 


from Eq. (12) to the case of three generations, we hnd agreement with the exact result: 


0.905 7-p|V.(M.)MMM.)W(M,)],| « 025GeV. 


(16) 


to be compared to the estimated value of 550 GeV. At hrst sight it may be surprising 
that the simple estimate outlined in Sec. leads to such a good agreement. However, 
given that there is quite a range possible for both Yukawa couplings and RH neutrino 
masses, it is in fact to be expected that, rather generically, one of the products of the 


12 








form (neutrino Yukawa coupling) x (RH neutrino mass) will in most cases dominate over 
the others, thereby mimicking the situation of only one RH neutrino being present. 

Glancing at Eq. (A-5b), one might be led to the assumption that, by raising the 
Majorana masses, the transition to negative values of can be pushed beyond Mgut, 
since below their mass thresholds the RH neutrino helds are integrated out. However 
such an attempt must fail since it is always overcompensated by the quadratic term in 
the RGE for m^, Eq. (A-5b), which is also illustrated by Fig. right panel. Only if we 
chose all Mj > Mqut could it be achieved, but at latest at that point some other new 
physics would probably appear which may completely change the situation. 

Having seen the examples in this section, it is evident that potentially non-trivial 
constraints can arise from avoiding the Z 2 breaking. However, depending on the context 
they may be stronger or weaker, so that the ultimate conclusion is that, for any given 
parameter point leading to Z 2 breaking, one must check whether or not the breaking may 
have a bad influence, in particular on the production of DM. 


5 Discussion & Outlook 

In this section, we will briefly discuss some further phenomenological implications and/or 
subtleties related to the parity problem. As we will see, none of the points changes the 
principal situation, however, the points discussed may lead to further constraints or they 
may be worth to be considered in a separate work. 


5.1 Unstable vacuum 


As we stated in Sec.|^ one can in principle obtain stronger constraints on the scalar sector 
if demanding that one of the Z 2 symmetric minima is also the global one. If this were not 
the case, the local minima would simply decay into the global vacuum. 

In repeating the analysis presented in Fig. (left panel) and checking for this additional 
feature, we obtain Fig. Here, the red points indicate that Z 2 breaking is the only 
constraint, while the blue points indicate that there is a deeper minimum emerging at a 
scale below the breaking scale A, thus giving even stronger constraints. As we can see, 
this is only relevant for light scalar masses. For large mass scales, Z 2 breaking evidently 
occurs before deeper minima can form. 

Unfortunately, it is rather difficult to infer any global qualitative statements from 


Eqs. (6a) to (6d) since they depend on the concrete values of the scalar couplings Ai, 2 , 3 , 4,5 
at a given scale. Intuitively, one would expect new minima to occur if m\ is either negative 
or at least small. For large physical scalar masses, m\ will be positive and large. Thus, 
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Figure 3: Re-run of the scan shown in Fig. (left panel), where blue points indicate that 
the Z 2 symmetric vacua are not the global ones. 


given that the scalar couplings have a tendency to run to smaller values for high scales, 
the scale of Z 2 breaking for growing approaches the scale where new minima emerge]^ 
Note that by a small mass parameter in this context we refer to the different combinations 
of mass parameters and quartic couplings, e.g. ^ in Eq. (6b). 


5.2 DM decay 


We have already mentioned that DM production may be affected by the parity problem. 
It is worth to note that the issue is however not so much about the mere abundance 
of the DM. While indeed a broken Z 2 symmetry would change the Feynman rules in the 
scotogenic model, and hence the amount of DM produced in the early Universe, this could 
typically be compensated by choosing slightly different parameter values which, given that 
neither the inert scalar mass nor the RH neutrino masses are known accurately, is not 
expected to be a major issue. 

The actual problem lies somewhere else: given that the radiatively induced VEV 
of r] is not tiny (i.e., it can be expected to be at least in the GeV to TeV ballpark for 
generic parameter choices) and that the coupling strengths are not extremely small, the 


®For example, the Higgs quartic coupling Ai is dominated by the top Yukawa coupling yt which gives 
a negative contribution, cf. Eq. (A-4a). 
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resulting DM candidate will generically decay very quickly. While decaying DM is not 
excluded by observations, its interactions typically need to be strongly suppressed e.g. by 
very heavy mediators or by tiny (~ gravitational) interaction strengths - see Ref. for 
a discussion of some example models. If this is not the case, the DM is in great danger 
to simply “decay away”, i.e., to have a lifetime smaller than the age of the Universe, so 
that it could not have survived until today when produced early in history. Potentially 
fast decays would furthermore tend to increase the freeze-out temperature, since decaying 
DM drops out of thermal equilibrium earlier than stable DM, so that on top of that (at 
least for non- relativistic DM particles) a strong modihcation of the abundance can be 
expected, too. 

One can easily see that the constraint arising from the lifetime is rather strong. For 
example, taking as Feynman rule for the vertex which allows an inert scalar to decay 
into two SM-like Higgses, the resulting decay width could be approximately computed 
to be F ~ v‘^\^/{IQurrir). If we conservatively estimate the DM lifetime to be at least 
larger than the age of the Universe (13.81 Gyr [^), corresponding to about 4.3 ■ 10^^ sec, 
we get an upper limit of about 1.5 ■ 10“^^ GeV on the decay rate. Even for a mass of 

= 500 GeV and a “small” VEV of = 1 GeV, the upper limit on the interactions 
strength would then be A < 6 ■ 10“^^. For a scalar potential like the one in Eq. ([^, this 
bound would extend to all three couplings A 3 , 4 , 5 , at least in the absence of very hne-tuned 
cancellations. Taking into account some information on the decay products, this limit 
would become even stronger [^. Similar estimates could be done in case one of the RH 
neutrinos was the DM particle. 

Hence, unless we resort to extremely small couplings or very strong hne tuning, a VEV 
of the inert scalar would quite generally destroy the ability of the scotogenic model to 
explain DM, simply because the interactions involved cannot be sufficiently suppressed. 


5.3 Thermal effects 


Thermal effects can be important for our analysis since broken symmetries will generally 
be restored at sufficiently high temperatures 36 . One might be led to the conclusion that 


this makes the entire discussion of Z 2 breaking superfluous, given that the breaking occurs 
at high scales only. This is, of course, too simple an argument since the real situation is 
much more complex. To understand this, let us reconsider the situation shown in Fig. [T] 
(right panel) with a 5 TeV heaviest RH neutrino and the lightest inert particle. We can 
capture the leading thermal effects by including the thermal corrections to the quadratic 
terms in the scalar potential 61 . This is achieved by a shift in the mass parameter 
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Figure 4: Effective thermal mass parameter for different temperatures T. 


m\ —)■ m\ + cT^, where T is the temperature and c is given by 69,72 : 


C — (s'! + ^gl) ^ 12 ’ 


(17) 


as can be calculated from the IR limit of the thermal scalar two-point function 75 77 


In this way, we obtain the effective thermal mass parameter which is shown in 

Fig. 1^ for different temperatures. Clearly, at high temperatures the Z 2 symmetry will 
be restored for a hxed renormalisation scale /r, as one would expect. However, since 
our DM candidate in this scenario is r]ji, which is a weakly interacting massive particle 
(WIMP) with a physical mass around 250 GeV, the temperatures relevant for the thermal 
decoupling will be much below M 3 , such that the thermal corrections have only little 
impact (cf. violet curves). The remaining question of how to choose the renormalisation 
scale can only be answered if one looks at the details of the DM production. We leave 
such details to be worked out in future studies. 

In combination with the discussion on DM decay, this shows that thermal effects will 
alter the details but not the overall picture: Z 2 breaking can occur and may have a 
signihcant impact on the phenomenology of the scotogenic model. Any phenomenological 
study should include this fact in their considerations. 
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6 Conclusions 


We have illustrated that the scotogenic neutrino mass model suffers from a parity problem, 
i.e., it is in danger that its intrinsic built-in Z 2 parity symmetry is broken by quantum 
effects driven by the heavy right-handed neutrinos. This could generate unwanted effects 
such as modifications of DM production. This issue imposes visible constraints on the 
parameter space available, in particular because the most generic solution, i.e., simply 
pushing the corresponding mass parameter in the Lagrangian to large enough values to 
avoid the breaking, does not work due to electroweak symmetry breaking being threat¬ 
ened. Thus, the scotogenic model suffers from tension from two different sides which 
considerably shrinks the allowed parameter ranges. 

After introducing the scotogenic model and its vacuum structure, we have presented an 
analytical approximation to compute the breaking scale, before exemplifying the resulting 
constraints numerically. Our considerations are based on the 1-loop renormalisation group 
equations of the scotogenic model which we have re-derived and, in passing, updated com¬ 
pared to previous versions in the literature. Summing up, we have revealed a somewhat 
subtle but non-trivial constraint on the scotogenic model which is able to strongly reduce 
the allowed parameter space. This makes it necessary for future considerations to check 
whether the parity problem exists for a certain choice of parameters, or not, to avoid the 
trap of studying physically irrelevant regions of the model. 
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A Renormalisation group equations 


The 1-loop RGBs for the scotogenic model have first been computed in Ref. 24 . We 


have re-derived those equations needed for the purpose of this paper, and have in passing 
taken the opportunity to update part of the earlier results. 

For convenience, we define the differential operator T) = The 1-loop RGBs 
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for the gauge couplings are those of a generic THDM 



= bigf (no sum!), 


(A-1) 


with b = (7, —3, —7). 

The quark sector of the scotogenic model is the same as that of the SM, such that the 
corresponding RGBs do not change]^ The RGBs for the leptonic Yukawa couplings are: 


= n + i/!*/! + T-^gf- , (A-2a) 

BA = ft I ?fttft + - HgA, (A-2b) 

where T^, = Tr (hl/i) and T = Ti + 3Y^Yu + . Bor the Majorana mass 

matrix, one hnds |^[7^ : 

T>M= |(hh^)M + M(hh^)*|. (A-3) 


Bor the quartic scalar couplings, we hnd the RGBs for a Z 2 symmetric THDM |^: 


VXi — 12 X 1 + 4A3 + 4A3A4 + 2A4 + 2A5 + - (gi + 2 glgl + Sg^) 

~ 3Ai (^gi + 3g2) + 4AiT — 4 T 4 , 

VX2 = I2A2 + 4A3 + 4A3A4 + 2A4 + 2A5 + - [gi + 2 ,g\g\ + Sg^) 

— 3 A 2 [gl + 3 g 2 ) + 4A2T!2 — 4T4,g, 

= 2 (Ai + A2) (3A3 + A4) + 4A3 + 2A4 + 2A5 + - ((74 — 2 glgl + Sg ^) 

~ 3A3 (774 + 3772) + 2A3 (T + Ty ) — 4 Tye , 

DA4 = 2 (A4 + A2) A4 + 8A3A4 + 4A4 + 8A5 + 3 g \ g "^ 

— 3A4 (774 + 3772) + 2A4 (T + Ty ) + 4 Tj^e) 

DA5 = A5[2 (A4 + A2) + 8A3 + I2A4 — 3 (774 + 3772) + 2 (T + Ty )], 

®Note the implicit changes in gi^ 2 ^ though, by virtue of Eq. (|A-1[). 


(A-4a) 


(A-4b) 


(A-4c) 

(A-4d) 

(A-4e) 
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where we have used the abbreviations T 4 = Tr {Y^YfY^Yf. + YY^Y^Y^Y^ + SY^Y^Y^Yd^, 

T 4 ^ = Tr (h^hh^h) and T^e = Tr {h^hY^^Y^). 

The scalar mass parameters obey the following RGEs: 


Vml 


6Ximl + 2 ( 2 A 3 + A4) ml + 


2r - 5 (9? + 39?) 


(A-5a) 


Vml — QX 2 ml + 2 ( 2 A 3 + A4) m\ + ml 


2r. - 2 (9? + 39?) 


(A- 5 b) 


2 = 1 


where the last term in Eq. (A-5b) is characteristic for a scalar field coupled to Majorana 


fermions (see, e.g., Ref. 81 This term is the crucial point of our study. Conveniently, 


since it is nothing but a trace, it is invariant under the transformation that diagonalises M, 
such that we do not have to perform this diagonalisation explicitly. Furthermore, the de¬ 
coupling of the heavy Majorana helds must be carried out by hand, since we are working in 
a mass-independent renormalisation scheme. To this end, we match the neutrino Yukawa 
couplings and Majorana mass matrices in the basis where M is diagonal (7^. Apart from 


these two quantities, the only appearance of the Majorana masses is Eq. (A-5b), where it 
suffices to remove the corresponding contribution. 
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